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Abstract. We present a study of ISOGAL sources in the "intermediate" galactic bulge (|/| < 2°, |6| ~ l°-4°), 
observed by ISOCAM at 7 and 15 /jm. In combination with near-infrared (I, J, K s ) data of DENIS survey, 
complemented by 2MASS data, we discuss the nature of the ISOGAL sources, their luminosities, the interstellar 
extinction and the mass-loss rates. A large fraction of the 1464 detected sources at 15 fim are AGB stars above 
the RGB tip, a number of them show an excess in ([7]-[15])o and (K s -[15])o colours, characteristic of mass-loss. 
The latter, especially (K s -[15])o, provide estimates of the mass-loss rates and show their distribution in the range 
10" 8 to 10" 5 M©/yr. 
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1. Introduction 

ISOGAL is a detailed mid-infrared imaging survey of the 
inner Galaxy, tracing the Galactic structure and stellar 
populations (Omont et al. 2003). It combines 7 and 15 
Urn ISOCAM data with IJK S DENIS data (Epchtein et 
al. 1994, 1997; Simon et al. in preparation). The main 
goals of ISOGAL survey are to trace the large scale in- 
ner disk and bulge structure using primarily red giants 
old stars; to study the corresponding stellar populations; 
to trace the age and mass-loss of AGB stars; to deter- 
mine the number of (dusty) young stars and to map the 
star formation regions through diffuse ISM emission and 
extinction. Multicolour near- and mid-infrared data are 
essential to analyse these features. 

The scientific results of ISOGAL are summarized in 
Omont et al. (2003). The analysis of the first ISOGAL 
field is detailed in Perault et al. (1996). The discussion of 
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a small field in the inner bulge {I — 0°, b = +1°) (Omont 
et al. 1999) has confirmed and characterized the remark- 
able capability of ISOGAL to detect the mass-losing AGB 
stars. Glass et al. (1999) analysed the ISOGAL fields in 
Baade's Windows of low obscuration towards the bulge. 
Most of the detected objects towards Baade's Windows 
are late-type M stars, with a cut-off for those earlier than 
about M3-M4. In the same fields, Alard et al. (2001) have 
determined the general properties of a well defined sam- 
ple of semi-regular variables by obtaining the MACHO 
lightcurvcs in V and R for a large fraction (~ 300) 
of the ISOGAL objects in Baade's Windows. Recently, 
Schultheis k Glass (2001) and Glass & Schultheis (2002) 
have further investigated the luminous M-type giants in 
these Baade's Windows fields, in the near-infrared mainly 
taken from the DENIS and 2MASS (Skrutskie et al. 1995) 
surveys. 

In this paper, we report the study of nine ISOGAL 
fields in the intermediate galactic bulge (— 1.5°<^<+1.6°; 
-3.8°<6<-1.0°, b = +1.0° (Fig. 1); with a total area ~ 
0.29 deg 2 ) and more than 2300 detected sources. The 15 
iim and 7 /im ISOGAL point sources have been combined 
primarily with DENIS IJK S data in the ISOGAL-DENIS 
Point Source Catalogue which gives magnitudes I, J, K s , 
[7], [15] (Schuller et al. 2003), - and also with 2MASS (J, 
H, K s ) and MSX (mid-IR, Price et al. 2001) - to deter- 
mine their nature and the interstellar extinction. Analysis 
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Fig. 1. The rectangles show the limits of the ISOGAL 
fields studied in this paper. 

at five near- and mid- infrared wavelengths of the ISOGAL 
sources, together with the relatively low and constant ex- 
tinction, shows that the majority of the sources are red 
giants with luminosities just above or close to the RGB 
tip and with a large proportion of detectable mass-loss. 
The various colour-colour and magnitude-colour diagrams 
are discussed in the paper, together with values of the lu- 
minosities and mass-loss rates. 

The outline of the paper is as follows : In Sect. 2, we 
present the ISOGAL and DENIS observations, as well as 
complementary 2MASS (and MSX) data. Sect. 3 describes 
the cross-identification of the mid-IR ISOGAL sources 
between 7 & 15 //to bands and then with near- infrared 
sources. In Sect. 4, we discuss the determination of inter- 
stellar extinction in the line of sight of the bulge fields 
using the isochrone fitting. Sect. 5 presents the derivation 
of bolometric magnitude (Mboi) and luminosity for each 
star in the bulge fields. In Sect. 6, we discuss the nature 
of ISOGAL sources. Sect. 7 is devoted to the derivation of 
mass-loss rates based on (K s -[15])o colour. In Sect. 8, we 
consider a preliminary estimation of the total mass-loss 
rate in the intermediate bulge. 

2. ISOGAL observations and complementary data 

ISOGAL observations with ISOCAM (see Cesarsky et al. 
1996 for a general reference to ISOCAM operation and 
performances) at 6 pixel field of view of nine intermedi- 
ate bulge fields with good quality of 7 and 15 //to data are 
used in our study (see Fig. 1). These observations were 
performed during 1996-1998 at 15 //to (filter LW3 : 12- 
18 /zm or LW9 : 14-16 //to), and at 7 [im (filter LW2 : 
5.5-8.5 (j,m or LW6 : 7.0-8.5 fim) (Table 1). For all the 



observed ISOGAL fields, we have used the data of the 
ISOGAL Point Source Catalogue (PSC, Version 1) de- 
scribed in Schuller et al. (2003), reproduced, for our nine 
fields, in a shortened version in Table 2 available at CDS 1 . 
They were derived with an improved data reduction with 
respect to ISO Archive data, with the use of CIA ISOCAM 
software and of a special source extraction (Schuller et al. 
2003). From the comparison between the different obser- 
vations of the field at I = 0°, b = +1°, we find very small 
differences between the different filters : LW2 - LW6 = 
+0.05 mag & LW3 - LW9 = +0.1 mag. Hence, no correc- 
tions are applied in LW6 and LW9 magnitudes to jointly 
discuss them with the main LW2 and LW3 data. The to- 
tal number of sources in the nine fields are 2232 and 1464 
at 7 and 15 //to, respectively, with a grand total of 2354 
sources. The corresponding average source densities are 
~ 8xl0 3 deg -2 and ~ 5xl0 3 deg -2 at 7 and 15 //m, re- 
spectively (see detailed density for each field as shown 
in Table 1). The 7 and 15 //to magnitude limits of the 
ISOGAL PSC are given in Table 1 and vary from 8.6 to 9.5 
at 7 //to, and from 7.5 to 8.5 at 15 //to. Generally, these 
ISO cutoff magnitudes should roughly correspond to a de- 
tection completeness level close to ~ 50%, and depend on 
the source density, so that these magnitude limits vary by 
about 1 mag among the different fields. The histograms of 
the 7 and 15 //to total source counts derived from the 6" 
ISOGAL observations are displayed in Fig. 2. This figure 
also displays the data (available at IPAC, version 1.2) of 
the MSX sources detected in these fields (108 at 8 //m, all 
with [7]iso < 7.0 and 16 at 15 /ito, all with [15])iso < 4.5). 
They show the large difference in sensitivity and in the 
number of detected sources between the two surveys. We 
have checked the consistency of ISOGAL and MSX 15 //to 
magnitudes for the 13 strongest sources (with [15] < 3.5, 
in order to avoid the Malmquist bias). The average dif- 
ference is [15]iso _ [15]msx = -0.06+0.43, where the large 
dispersion is probably related to the source variability. 

In the near-infrared, the data mostly used in this paper 
are also those of the PSC of the ISOGAL-DENIS survey. 
They are complemented by the 2MASS data (from the 
Point Source Catalog of the 2MASS second release avail- 
able at IPAC and CDS). However, these 2MASS data do 
not fully cover all the nine ISOGAL fields. There are large 
blank areas (because of saturating stars) in three fields, so 
that only about 88% of the total area of the nine ISOGAL 
fields is covered. All sources in these 2MASS data have 
detections in the three JHK S bands. The DENIS near- 
infrared data are acquired from the DENIS survey with 
special observations of the Galactic bulge (Simon et al., 
in preparation) in the three bands, I (0.80 //to), J (1.25 
//to) and K s (2.17 //to), which completely cover the nine 
considered ISOGAL bulge fields. Therefore, we have taken 



1 Table 2 contains the catalogue of ISOGAL-DENIS sources 
from nine bulge fields with their 2MASS & MSX counter- 
parts. The catalogue will be available via the VizieR Service 
at the Centre de Donnees Astronomiques de Strasbourg (CDS, 
http://vizier.u-strasbg.fr/viz-bin/VizieRl. 
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Fig. 2. 7 and 15 fim source distributions in half magnitude bins. ISOGAL magnitudes are converted into flux densities 
using the relations given in Schuller et al. (2003). The dotted lines show the MSX A (~ 8 fim) and D (<~ 15 fim) 
sources cross-identified with ISOGAL 7 and 15 fim sources. MSX magnitudes are derived using the zero magnitude 
flux given by Egan et al. (1999). 



the more complete DENIS data as the main reference for 
near-infrared associations. 

The histograms of all the DENIS I, J, K s sources in 
the observed fields are shown in Fig. 3. The complete- 
ness limit is close to 12.0 in K s band, 14.0 in J band and 
15.5 in I band, respectively. However, only near-infrared 
sources with K s & 11.0 are used for near- infrared cross- 
identifications, in order to limit the number of spurious 
associations (see Sect. 3 and Table 1). Fig. 3 also dis- 
plays the histograms of J, H and K s 2MASS sources. We 
have compared the DENIS and 2MASS photometry of 
1963 stars common in ISOGAL-DENIS PSC and 2MASS 
in our nine fields. We find small differences in the aver- 
ages: K D enis-K 2 mass = -0.00±0.07 and Jdenis~J2MAss 
= -0.02±0.09; hence no corrections are applied. 

Table 1 shows the used ISOGAL, DENIS and 2MASS 
observations in detail. It should be noted that the major- 
ity of the observations were performed on dates with large 
differences (from two months to three years), except for 
7 and 15 fim observations performed on the same date in 
six fields out of nine, and the DENIS, 2MASS observa- 
tions performed on the same date for all sources in their 
respective bands. Such differences in dates are comparable 
to the period of Long Period Variables (LP Vs) . One thus 
expects that the corresponding colours of the few large 
amplitud & variables (Miras) may be strongly affected. 
This uncertainty will propagate through to the final un- 
certainties in luminosity and mass-loss as well. From the 
detailed studies of LPVs in the two Baade's Window fields 



(see Glass et al. 1999, Alard et al. 2001, Schultheis & Glass 
2001, Glass & Schultheis 2002), one can estimate that the 
variability r.m.s. error in the near-infrared is ~ 0.3 mag 
for strong variables (miras), i.e. most of the strong sources 
with [15] < 4, and less than 0.1 mag for most others. In 
order to reduce this uncertainty by about a factor \/2, we 
have systematically used: 1) The average of DENIS and 
2MASS J and K s magnitudes when available; however, 
the DENIS sources are saturated at K s <i 6.5, and, for 
such few strong sources, only 2MASS data are used (Sect. 
3). 2) The average of ISOGAL and MSX 15 fim magni- 
tudes for the few (13) strong sources detected by MSX 
with [15]iso < 3.5. 

3. Cross-identification of ISOGAL and 
near-infrared sources 

In order to perform the association between 7 and 15 fim 
sources in the ISOGAL PSC, first a global offset between 
the two sets of coordinates has been determined, then dis- 
tortion coefficients have been computed to correct for the 
rotation effects (Schuller et al. 2003). Finally, 7 and 15 fim 
associations were searched within a radius of twice the size 
of the pixel (2x6"= 12"). Only - 8% (122) of the 15 fim 
sources have no 7 fim counterpart. However, out of these 
122 15 fim sources, a large fraction (77% : 94 sources) have 
been associated with K s (K s <> 11.0) DENIS or 2MASS 
sources, which reflects the incompleteness of weak 7 fim 
sources. The r.m.s. of the [7]-[15] separations has a small 
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Table 1. Log of ISOCAM, DENIS and 2MASS observations in the intermediate bulge fields 
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value equal to 1.5", which shows the good quality of most 
of the associations. Indeed, 98% have a quality factor 3 or 4 
in the ISOGAL PSC which warrants almost certainly real 
association (Schuller et al. 2003), only 2% have a quality 
factor 2 which means that the association is probably real, 
but should be considered with some care (a few doubtful 
associations with quality factor 1 have been discarded). 

A cross-correlation algorithm similar to the one used 
to associate 7 and 15 fim data has been used to look for 
DENIS counterparts to ISOGAL sources (Schuller et al. 
2003). The density of DENIS sources is so high that the 
DENIS catalogues were cut to 36000 sources per square 
degree (K s ^ 11.0, see Table 1). A search radius of 7" was 
initially used for the cross-correlation between ISOGAL 
and DENIS sources, and the association was limited to 
DENIS sources detected in K s . Figure 4 shows the his- 
togram of positional difference of ISOGAL-DENIS cross- 
identified sources in our nine fields. The r.m.s. of the dis- 
tances of ISOGAL-DENIS association has a value equal 
to 1.4". About 98% of the ISOGAL sources (either 7 or 
15 fj,m) have thus initially an association with a DENIS 
source with K s <J 11.0. Even with such a large search ra- 
dius, because the proportion of real associations is very 
large and only associations with the smallest separation 
are retained, the fraction of false associations is expected 
to remain ~1%. In addition, the quality of the associations 
is characterised by a specific flag and we have dropped in 
Table 2 the associations with the poorest quality (see be- 
low). Such a large radius was indeed needed since, in some 
cases, real associations have such large separation, because 
of the large pixels and undersampling of ISOGAL, in par- 
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Fig. 3. DENIS source counts in I, J & K s bands in half 
magnitude bins (thick continuous lines). Thin continuous 
lines show the DENIS sources detected in J and K s bands. 
Dotted lines show the 2MASS sources detected in J, H & 
K s bands. The used 2MASS data (2nd release) contain 
only sources detected in the three JHK S bands and cover 
only ~ 88% of the total area of the considered ISOGAL 
fields. Long-dashed line shows the average K s magnitude 
cut of DENIS sources (K s < 11, see Table 1 and Sect. 3) 
cross-identified with ISOGAL sources. 
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Fig. 4. Histogram of positional difference (in arcsec) be- 
tween ISOGAL and DENIS cross-identified sources. Solid 
line shows all ISOGAL/DENIS sources. Dotted line shows 
the sources having ISOGAL/DENIS association quality 
flags 4 & 5. Dashed line shows the sources having quality 
flags 2 & 3. Dotted-dashed line shows the sources having 
quality flag 1. Associations with quality flag 1 and those 
with a distance larger than 5"have been finally discarded. 

ticular in the cases of blended or slightly extended sources. 
This is extensively discussed in Schuller et al. (2003). 
About 1% of the ISOGAL sources (23) have a DENIS 
counterpart which are saturated (K s ^ 6.5, see Fig. 7), out 
of which 18 saturated DENIS sources have been replaced 
by the 2MASS sources (as indicated in Table 2 available at 
CDS). The five other sources are not identified in the pub- 
lished 2MASS fields. Out of these 18 sources, 50% sources 
are identified as foreground and taken out of the stars 
further discussed in the paper (see text, Sect. 4) and the 
remaining are shown in various figures (see Figs. 8, 9, 11 
& 12). Most of the ISOGAL sources cross- identified with 
DENIS ones (92%) have a good quality association (qual- 
ity flags 4 & 5; see Schuller et al. 2003 for details), 7% 
sources have more doubtful association (quality flags 2 & 
3) and only less than 1% ones have very bad association 
(quality flag 1). For the present study, we have dropped in 
Table 2 the sources which have association quality flag of 
1 (15 associations) and those with associated distance of 
more than 5" (12 associations, mostly with quality flag 2). 
The number of contaminating false matches is thus kept 
well below 1%. The respective numbers of 7 fim alone, 15 
fim alone & 7+15 \im sources not associated at all with 
DENIS sources are 21, 32 & 18. However, most of these 
7+15 \im sources are not in the DENIS catalogue because 
they are too strong and above the saturation limit; 13 of 
them are present in the 2MASS data. Most (18) of the 
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Fig. 5. Colour-magnitude diagrams (J-K s /K s ) of DENIS 
sources in the ISOGAL intermediate bulge fields. The four 
isochrones (Bertelli et al. 1994), placed at 8 kpc distance 
for a 10 Gyr population with Z=0.02 are shown in each 
figure for Ay = mag, for the Ay limit adopted for the 
"foreground" sources (defined as "F" in top left box), for 
mean Ay values of the field and for the Ay limit adopted 
for the "background" sources (defined as "B"), respec- 
tively. 

sources of 7 /jto alone are also associated with 2MASS 
sources with K s ^ 11.0, so that the proportion of 7 [im 
sources not associated with a near-IR source is negligible 
(~ 0.4%). On the other hand, only four sources of 15 \im 
alone are associated with 2MASS sources with K s ^ 11.0; 
it results in a more important, but still small, fraction (~ 
2%) of 15 /zm sources without near-IR associations, espe- 
cially for the sources detected only at 15 (im (~ 23%). The 
latter could be young stellar objects (YSOs, Sect. 6.2); 
however, half of them have poor quality flags which could 
be characteristic of the spatial extension of YSOs (Schuller 
et al. in preparation) or cast a doubt on the reality of unas- 
sociated sources. Finally, taking into account the 15+12 
rejected associations, the percentage of ISOGAL sources 
with a near-infrared association is larger than 97%. 

4. Interstellar extinction and foreground disk stars 

4.1. Interstellar extinction 

The adjunction of DENIS and 2MASS near-infrared data 
adds much to ISOGAL data, by providing different and 
more sensitive colour indices, as well as estimates of the 
interstellar reddening. The K s /J-K s magnitude-colour di- 
agrams of DENIS sources in the bulge fields show a well- 
defined red giant sequence shifted by fairly uniform ex- 
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tinction (see Fig. 5), with respect to the reference K s o vs 
(J-K s ) of Bcrtelli et al (1994) with Z=0.02 and a dis- 
tance modulus of 14.5 (distance to the Galactic Center : 
8 kpc). The near- infrared colours of this isochrone have 
been computed with an empirical T cf f-(J-K) colour rela- 
tion built by making a fit through measurements for cool 
giants (see Schultheis et al. 1998; 1999). We have assumed 
that Aj/Ay = 0.256; A Ks /A v = 0.089 (Glass 1999). The 
average value of the interstellar extinction (Ay) for each 
bulge field is shown in Table 1 and in Fig. 5. The mean 
Ay value for each field has been determined by a gaussian 
fit to the Ay distribution. 

The DENIS sources with anomalously low values of 
Ay are probably foreground. They are visible in Fig. 5 in 
each of the bulge fields (stars around the isochrone with 
Ay <~ mag and stars clearly left of the bulk of the stars 
grouped around the isochrone with mean Ay of the field) . 
For each field we define an isochrone "F" (Fig. 5) for which 
we assume that all the sources left of it are foreground 
stars which will be no longer considered in the following 
discussions of "bulge" stars. 

There are also a number of stars with J-K s values sig- 
nificantly larger than the bulk of the other stars in each 
bulge field (right of the isochrone defined as "B" in Fig. 5). 
We can see three reasons for such an excess in J-K s : 1) in- 
trinsic (J-K s ) excess induced by a large mass-loss, which 
should be accompanied by a large 15 /zm excess; 2) excess 
in Ay which should probably be due to a patchy extinction 
on the bulge line of sight (additional extinction from dust 
layers behind the Galactic Center for background stars ap- 
pears unlikely at such high galactic latitudes); 3) spurious 
association or wrong photometry which is rather unlikely 
for DENIS sources well above the detection limit. 

For the stars with such large values of J-K s , right of 
the isochrone "B" in Fig. 5, we distinguish: a) bright stars, 
with K s < 8, which can be bulge AGB stars with large 
mass-loss, for which we again use the mean extinction cor- 
rection Ay of the field; and b) faint stars, with K s > 8, for 
which we assume low mass-loss and large extinction and 
we determine their specific extinction from J-K s and the 
zero-extinction isochrone as for the bulk of bulge sources. 

To summarize, for the following discussions where 
dereddening is essential - colour-magnitude diagrams, lu- 
minosities, mass-loss rates, we have applied the following 
prescriptions : We have discarded all sources left of curves 
"F" in Fig. 5. For the others, with J and K detections, 
besides the few exceptions mentioned above, we have de- 
termined their specific extinction from J-K s [with Glass 
(1999) values for Aj/Ay and Ak s /Av] and the quoted 
zero-extinction isochrone of M giants. For all the others, 
especially with no JK associations, we have used the mean 
extinction Ay of the field (Table 1). The extinction ratios 
A[7]/Ay and A[ 15 ]/Ay in the ISOGAL bands are still un- 
certain (Hcnnebelle et al. 2001, Jiang et al. 2003). We have 
used the values A [7] /A v = 0.020; A [15 ]/A v = 0.025 rec- 
ommended by Hennebelle et al. for "disk" fields; however, 
the extinction corrections remain small with the low val- 



ues of Ay considered, and the effects of the uncertainty 
on A[7] and A[ 15 ] are practically negligible. 

4.2. Proportion of foreground and background disk 
stars 

It is clear that the method discussed above for identifying 
foreground stars from the value of J-K s and the inferred 
extinction, fails to recognize foreground disk stars located 
between the last dust layers and the bulge. Similarly, the 
extinction gives no way of distinguishing background disk 
stars, beyond the bulge, from bulge stars. 

In order to derive an order of magnitude of the propor- 
tion of such foreground and background disk stars, mixed 
with bulge stars with similar extinction, one can use any 
reasonable simple Galactic model for the M-K giants of 
the bulge and the disk, e.g. the one by Wainscoat et al. 
(1992). In Appendix A, we have thus estimated the ap- 
proximate fraction of counts expected from the bulge, the 
disk r < Ro, and the disk r > R , as being about 65%, 
23% and 12%, respectively, on a typical line of sight for 
our fields, I — 0, b — 2° (R is the distance to the Galactic 
Center) . 

Indeed, stars of the actual bulge and of the central 
disk arc indistinguishable in our analysis. We should thus 
consider as really foreground and background, only disk 
stars outside of the bulge radius Ri = 2 kpc of Wainscoat 
et al (1992). Their fractions are only about 14% and 6% 
respectively. However, part of the foreground stars have 
already been identified from the low value of their extinc- 
tion (Sect. 4.1) and taken out of the stars further discussed 
below. Their number is 110, i.e. about 8% of the total. 

To summarize, one may expect that the sample of 
ISOGAL sources further discussed below contains about 
6% of foreground disk stars, 6% of background disk stars, 
and 88% of central stars within the bulge limits - ~70% 
of bulge stars and ~18% of central disk stars. However, 
these percentages must be considered as indicative since 
they vary with the Galactic coordinates and they result 
from approximate estimates with a sketchy axisymmetri- 
cal bulge model. 

5. Luminosities 

After dereddening (Sect. 4.1), the bolometric magnitudes 
(Mboi) are derived for the stars using the bolometric 
corrections (BCk) and a distance modulus of 14.5. We 
have used an approximate relation (Fig. 6) between BCk 
vs (K s -[15])o to derive the bolometric corrections for 
the ISOGAL sources by combining the BCk values de- 
rived as follows from Frogel & Whitford (1987) and from 
Groenewegen (1997). For (K s -[15])o < 3, we have used 24 
ISOGAL Baade's Windows sources (shown by open cir- 
cles in Fig. 6) which are common in Frogel & Whitford's 
(1987) catalogue with a value for BCk, to derive the av- 
erage curve displayed in Fig. 6. We have extended this 
curve for (K s -[15])o > 3 from the approximate relation 
between BCk vs (K-[12]) for IRAS AGB stars proposed 
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by Groenewegen (1997). For this purpose, we have used 
an approximate empirical relation between y = (K-[12]) 
and x = (K s -[15])o - namely, y = x - (x/7) 3 - derived 
from the modelling of infrared emission of nearby oxygen- 
rich AGB stars (Groenewegen, private communication). 
The BCk values from Groenewegen (1997) are shown as 
open squares in Fig. 6. However, they have been slightly 
rescaled by adding 0.25 mag to match them with Frogel & 
Whitford's values for (K s -[15])o 5s 3. Finally, a polynomial 
fit has been performed between BCk and x = (K s -[15])o 
(see caption of Fig. 6). We have checked that such a sim- 
ple approximate determination of Mboi agrees within a few 
tenths of magnitude with an elaborate integration of all 
DENIS and ISOGAL flux densities along the method of 
Loup et al. (in preparation). 

The luminosity of each star is derived from Mboi us- 
ing the classical relation, L(L©) = i()-( M boi-4.75)/2.5 5 and 
is available in Table 2. Figs. 7a, b show the histograms of 
bolometric magnitudes and luminosities of the sources in 
ISOGAL bulge fields. These distributions are clearly in- 
complete below the luminosity of the bulge RGB tip. One 
should keep in mind that for ~ 12% of the sources, the lu- 
minosities thus derived are wrong because their distances 
are significantly different from the bulge distance, either 
foreground or background on the line of sight. This will 
spuriously broaden the luminosity distribution. In partic- 
ular, most of the sources appearing with very large lumi- 
nosities are probably foreground stars with smaller lumi- 
nosities. In addition, even for the other sources within the 
"bulge" limits (6-10 kpc), the difference in distances with 
respect to 8 kpc induces a spread by about one magnitude 
(but with an r.m.s. <, 0.2 mag because the distribution 
is very centrally peaked). Other uncertainty terms, to be 
combined in quadrature, include the residuals to the BCk 
fit and the BCk calibration (r.m.s. 5s 0.2 mag), uncer- 
tainties in the K and 15 fim photometry (r.m.s. 5s 0.2 
mag), dcrcdddcning for large extinctions (r.m.s. < 0.2 
mag), variability (r.m.s. <~ 0.2 mag for strong variables, 
< 0.1 mag for others). The combined r.m.s. uncertainty 
for the luminosity of bulge sources is thus probably less 
than 0.4-0.5 mag. 



6. Nature of the ISOGAL sources 

Figs. 8, 9 & 11 show the (K s -[15]) /[15] , (K s -[15]) /K s0 
and [7]- [15]/ [15] colour-magnitude diagrams of ISOGAL 
sources, respectively. The first two diagrams are more di- 
rectly related to the luminosity and mass-loss derivations 
(Sects. 5 & 7), and luminosity and mass- loss rates estimate 
are displayed in Figs. 8 & 9. However, such derivations 
rely on a good extinction correction of K s magnitudes. 
On the other hand, the diagram [7]-[15]/[15], which de- 
pends little on extinction, provides a straightforward view 
of the nature and the amount of dusty sources of various 
classes. It even allows a direct measure of the mass-loss 
rates, but with less sensitivity, especially for the largest 
and the smallest ones. 




2 4 

(K s -[15]) 

Fig. 6. Approximate relation between bolometric correc- 
tion, BC K , versus x = (K s -[15]) (BC K =2.982+0.0962x- 
0.1653x 2 + 0.1499x 3 -0.04173x 4 +0.003351x 5 ) derived us- 
ing the BCk values from Frogel & Whitford (1987) (open 
circles) and Groenewegen (1997) (open squares) (see in 
the text). 
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Fig. 7. Histograms showing bolometric magnitudes and 
luminosities of the sources in ISOGAL bulge fields. Solid 
lines show all ISOGAL sources detected at K s and 15 \im. 
Dotted lines show the sources having ISOGAL/DENIS as- 
sociation quality flags 4 & 5. Short-dashed lines show the 
sources having quality flags 2 & 3. Long-dashed vertical 
lines indicate the bulge RGB tip (M bo i <~ -3.5; L ~ 2000 
L ). 



Figs. 8, 9 & 11 are very similar to the same diagrams 
in Omont et al. (1999) for the field at I = 0, b = +1° (see 
also Glass et al. 1999). One can distinguish three main 
classes of sources, corresponding to three regions delimited 
in Figs. 8, 9 & 11 : 

There are first, a number of stars with luminosity 
(< 2000 L ) and colours [(K s -[15]) < 1] and [7]-[15] con- 
sistent with dustless stars on the RGB close to its tip, 
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(K s -[15]) 

Fig. 8. [15]/K S -[15] magnitude-colour diagram of 
ISOGAL sources with DENIS/2MASS counterparts. 
The open circles show the sources having ISOGAL- 
DENIS association quality flags of 4 & 5. The crosses 
show the sources having association quality flags of 2 & 
3 (see the text). The filled hexagons display the ISOGAL 
sources associated with 2MASS counterparts (since the 
DENIS sources [K s < 6.5] are saturated, they have 
been replaced by the 2MASS sources). The remaining 
saturated DENIS sources (which have no published 
2MASS counterparts) are shown by filled squares. The 
bright ISOGAL sources associated with 2MASS sources 
(which have no DENIS counterparts) are shown by open 
triangles. Sources without 7 \xm, associations are displayed 
by filled circles. The filled triangles show the ISOGAL 
sources associated with 2MASS sources but having no 
DENIS counterparts. The IRAS sources are shown by 
open squares. However the four very red IRAS sources 
11-4 (Table 3) are outside the limits of the figure. The 
asterisks show the sources detected at 8 [ira in the MSX 
survey which are associated with strong ISOGAL sources 
([15] < 4.5 mag). Among them the dark ones are those 
associated with the MSX D (15 fim) sources ([15] < 
3.5 mag). The curves from bottom to top represent the 
luminosity estimates from 500 to 20000 L , respectively 
(see Fig. 5). Various regions discussed in the text (see 
Sect. 6) are indicated with boundaries as dotted lines. 
The approximate scale of mass-loss rate displayed at the 
top of the figure is derived from Eq. (1) given by Jeong 
et al. 2002 (see Sect. 7). 

which is at Ko ~ 8.0 as quoted by Frogel et al. (1999) 2 , 
which corresponds to Mboi — -3.5 and L ~ 2000 L Q . Their 

2 As quoted in footnote 2 of Omont et al. (1999), we neglect 
the small difference of K-K s ~ 0.04, for such stars. 



spectral types should be ~M3-M5 (Frogel & Whitford 
1987, Glass et al. 1999), which corresponds to colour 
(K s -[15]) ^ 0.71-0.78 with T cff ~ 3600-3400 K, respec- 
tively. 

Then, a whole sequence of sources display luminosities 
slightly above the RGB tip (2000-5000 L Q ) and colours 
with a weak 15 \xm excess [1 < (K s -[15])o <2], character- 
istic of small amounts of dust and thus of weak mass-loss; 
they were called "intermediate AGB" stars by Omont et 
al. (1999). The K s0 magnitude range from the base to tip 
of the sequence (Fig. 9) corresponds to M spectral types 
from about M6 to M8-M9 (Frogel & Whitford 1987, Glass 
et al. 1999). It is known from Alard ct al. (2001) that 
they are semi-regular variables with relatively short peri- 
ods and weak amplitudes. 

Finally, less numerous sources, in the top right cor- 
ner of Figs. 8 and 11 ([15] < 4.5, (K s -[15]) > 2, 
([7]-[15])o 1.3), may have larger luminosities and super- 
wind mass-loss characteristic of Miras (see Glass et al. 
1999, Alard et al. 2001, Omont ct al. 1999). Most of them 
are detected at 8 \xm in the MSX survey (Price et al. 2001) 
as seen in Table 2. 

6.1. Brightest and Reddest ISOGAL sources 

We have identified 16 IRAS sources within the fields of 
Table 1. They are associated with 15 bright ISOGAL stars 
([15]o ^ 6.0). Table 3 shows the brightest and reddest 
ISOGAL sources with IRAS and MSX associations when 
available. These sources are shown by open squares in 
Figs. 8, 9, 11 & 12. Out of these 15 sources, the four bright- 
est ones ([15]o < 1.5) are denoted by "11-4" in Table 3 and 
Figs. 9 & 12. These four sources visible in Figs. 9 and 11 
are not shown in Figs. 8 & 12. All of them have extremely 
red values of (K s -[15])o characteristic of large amount of 
dust. Their [12]-[25] IRAS colours are characteristic of 
AGB stars and not of young stellar objects or post- AGB 
stars. They are thus probably AGB stars with very large 
mass-loss rates, > 10 -5 M Q /yr (see Sect. 7 and Table 3). 
Table 3 displays approximate values of their bolometric 
fluxes computed from an integration of their available flux 
densities, either following the prescriptions of Loup et al. 
(in preparation) or from a more approximate integration 
when (K s -[15))o is too large to use Loup et al's method. 
All of them correspond to large values of bolometric lu- 
minosities if the sources are at the Galactic Center dis- 
tance. Such luminosity values are consistent with the very 
large mass-loss rates and both quantities point out to rel- 
atively large values of the initial mass. However, it is not 
clear whether these stars belong to the central bulge/disk 
or to the foreground disk (the four of them have \b\ ~ 
1). Their large intrinsic colours make it very difficult to 
determine their interstellar extinction. Their following in- 
dividual properties may help to decide about their bulge 
membership: 

- II - It is the brightest 15 \im source and its ISOGAL 
magnitudes are uncertain because of saturation. It be- 
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Table 3. Catalog of brightest and reddest ISOGAL sources. Data for each source are displayed in two lines, with 
the ISOGAL standard name (e.g. ISOGAL-PJ173812. 5-293938) and the IRAS name if any. In order to make the easy 
comparison with ISOGAL 7 & 15 tim magnitudes, [7] & [15], and IRAS flux densities (in Jy) at 12 & 25 fim, S12 & 
S25, MSX intensities are given in magnitudes for bands A (<~ 8 /xm) and D (~ 15 urn), and in flux densities (in Jy) 
for bands C (~ 12 iim) and E (~ 21 iirri). 



Name (ISOGAL-P..) 




J 


K s 


[7] 


[15] 


A v 


Vlsr 


S12 


S25 


T 2 


M 


Id. 


IRAS Name 


b 






msxA 


msxD 






msxC 


msxE 










(deg) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(km/s) 


(Jy) 


(Jy) 


10 L 


(M Q /yr) 




J173812. 5-293938 


-1.5 


11.75 


8.16 


5.26 


2.89 


6.3 








0.32 


4.9xl0~ b 


AY 


J174122. 7-283146 


1.0 
-0.1 


>14.16 


11.80 


3.47 


1.54 


8.9 


-51.0 


2.0 


4.3 


0.89 


3.2 xl0~ 5 


13 


17382-2830 


1.0 






3.69 


1.87 






2.5 


4.7 








T1 7/) /(CO C 0*71 /I /I O 

J 1 / 4400.0-2 ( 1442 


1 A 

t.4 


1 o oo 
12. 06 


i n on 
1U.80 


U.90 


-0.51 


6.1 


-12. o 


i k n 
lo.U 


51.3 


a.oi 


4.1 X 1U 


en 


17418-2713 


1.0 






1.21 


-1.28 






31.2 


80.3 








J174506.1-271516 


1.4 


10.92 


8.67 


8.43 


5.84 


5.4 








0.16 


2.9xl0~ 7 


Y 


J174522. 9-271109 


1.0 
1.5 


>15.88 


12.50 


2.65 


1.42 


8.0 




4.5 


4.9 


3.67 


3.6 xl0~ 5 


14 


17422-2709 


1.0 






2.74 


1.23 






4.6 


5.7 








(J174533.2-304943) 


1.6 


(11.13) 


(7.64) 


(6.66) 


(6.13) 


(11.9) 






1.6: 


0.66 




YY 


17423-3048 


-1.0 






(6.07) 


















J174548. 6-304312 


-1.5 


>15.26 


10.14 


4.38 


1.58 


12.3 




3.4 


5.3 


0.17 


1.7xl0 -6 


12 


17426-3042 


-1.0 






3.67 


1.85 






2.7 


4.2 








J174939.4-292723 


0.0 


12.31 


9.85 


6.34 


4.16 


7.6 






2.5 


0.12 


8.1xl0~ 6 


15 


17464-2926 


-1.0 

























1 value comes from OH maser detection; 

2 if one assumes a distance of 8 kpc; 

3 the ISOGAL magnitudes and M bo i of II are uncertain because the ISOGAL flux densities are 
above saturation limits of the detectors. 



longs to the most extreme class of OH/IR stars with a 
massive cold circumstellar envelope as shown by its ex- 
ceptional [12]-[25] colour and strong 10 and 18 [im ab- 
sorption bands in the IRAS LRS spectrum (IRAS Science 
Team 1986). The radial velocity and the very large bolo- 
metric flux point to a foreground disk object. 

12 The bolomctric flux favours bulge membership. 

- 13 - It is also an OH maser, but its radial velocity 
and relatively small bolometric flux point to bulge mem- 
bership. 

- 14 - The bolometric flux and the blue [7]-[15] colour 
favour foreground disk membership. 

Three other stars with very red [7]-[15] colours (Fig. 
11) are also listed in Table 3. Their nature - AGB, post- 
AGB, PN, YSO - is uncertain, and they are not detected 
by MSX. We have also added an IRAS source, YY, with 
no ISOGAL point source counterpart: 

- 15 - From its position in the [7]/[7]-[15] diagram, it 
could be either a young stellar object (YSO) or an AGB 
star (Schullcr et al. in preparation). However, its non de- 
tection at 12 [im together with the large value of the ratio 
S25^m/Si5 M m ~ 6, and its small bolometric flux could in- 
dicate a young stellar object (or a planetary nebula). 

- AY - Again, from its position in the [7]/[7]-[15] and 
K s0 /(K s -[15])o diagrams, and from its bolometric flux, it 
could be a YSO or an AGB star. 

Y - It could appear as a good YSO candidate. 
However, its ISOGAL magnitudes could have been af- 



fected by the track of the saturated star II in the ISOGAL 
images. Its completely anomalous colour (K s -[7])o = -0.15 
also points to a wrong 7 fim magnitude. 

YY This source, IRAS17423-3048, is probably 
a massive YSO from its IRAS 25-60 fim flux densities. 
It has no convincing association in the ISOGAL-DENIS 
PSC, although the ISOGAL-DENIS-MSX source dis- 
played within brackctts in Table 3 is within the limits of 
a possible association. On the other hand, there is some 
extended emission at its position in the 15 lira ISOGAL 
image. 

6.2. Other detected sources 

All other IRAS detected sources have [12]-[25] colours 
indicative of AGB stars. The fact that there are prac- 
tically no other stars in Fig. 11 meeting the criteria 
([7]-[15] > 1.8 with [15]<4.5) for credible young star can- 
didates defined by Felli et al. (2000, 2002), is not surprising 
in such high galactic latitude lines of sight. 

As quoted in Glass (1999), much of the scatter 
(~ 1 mag) seen in the "AGB sequence" of Figs. 8 & 11) 
may be due to the distribution in the depth of the bulge 
and the central disk (Glass et al. 1995). Furthermore, in 
Figs. 8 and 9, it seems unlikely that there is in addition 
an anomalously large number of stars far above this se- 
quence by up to ~ 1.5 mag. Such stars with anomalously 
bright ISOGAL magnitudes with respect to their colours, 
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appear more clearly in the [7]-[15]/[15] colour/magnitude 
diagram of Fig. 11. It seems that most of them could not 
be bulge AGB stars with very high luminosity, because of 
the age of the bulge population and of the low mass-loss 
rates deduced for most of them from (K s -[15]) or [7]-[15]. 
Most of them are probably foreground AGB stars with 
moderate luminosity and mass-loss. Their population and 
the distribution of the values of the 15 [im magnitude 
excess with respect to the bulge sequence look in rough 
agreement with the proportion and the distance distri- 
bution expected for the disk foreground stars (Sect. 4). 
Similarly, most sources which lie below the main bulge 
AGB sequence in Figs. 8, 9 & 11, are likely background 
disk AGB stars. 

There are a few stars (Fig. 9) with large values of 
K s -[15]) (> 1) below the RGB tip. But their case is not 
clear. Only a few of them have values of ([7]-[15])o > 0.6. 
Most of the sources have poor photometry because they 
are close to the 15 [ira detection limit or of poor asso- 
ciation quality. In the following discussion, we will drop 
the most doubtful cases, i.e. all stars with K s o > 9.0 (Fig. 
9). The nature of the few remaining cases is unclear. A 
number could be background stars with large mass-loss. 
Additional observations are required to confirm the nature 
of these sources and to exclude the possibility of spurious 
associations with the K s sources. However, such a 15 /im 
excess could be an indication that some RGB stars close 
to the RGB tip have a significant mass-loss, as found in 
globular clusters by Origlia et al. (2002). 

The 28 sources detected only at 15 [im (Sect. 3) 
are absent from all colour-magnitude diagrams. However, 
they remain candidates for YSOs or Planetary Nebulae, 
especially for half of them with relatively good qual- 
ity flags. The strongest one, ISOGAL-PJ1753178-280433, 
[15] = 5.88, is identified in SIMBAD as a Planetary 
Nebula JaSt 90 (Van de Steene & Jacoby 2001). 

Fig. 12 shows [7]/K s -[7] magnitude-colour diagram of 
ISOGAL sources. The [7]/K s -[7] diagram seems to be the 
best criterion for the detection of large amplitude LPVs 
(Glass et al. 1999, Schultheis et al. 2001). The region of 
high mass- loss AGB stars ((K s -[7]) > 1, [7] < 7) is ex- 
plicitly drawn in this diagram. 

7. Mass-loss rate (M) determination 

It is generally agreed that the colour (K-[12]) (where [12] 
is the magnitude corresponding to the IRAS 12 [im flux 
density) is a good indicator of the mass-loss rate M of AGB 
stars (see e.g. Whitelock et al. 1994, Le Bertre & Winters 
1998, Jeong et al. 2002). It is clear that (K s -[15]) should 
be an equally good gauge of M, with a slightly different cal- 
ibration. However, the relation between M and (K s -[15])o 
or (K-[12]) directly depends on the nature of dust and on 
the gas-to-dust ratio. For instance, there is evidence from 
the preliminary analysis of the ISOGAL CVF data of a 
few of these AGB bulge sources that their dust is domi- 
nated by alumina for low and intermediate mass-loss rates 
(Blommaert et al. 2001 and in preparation). In addition 



Log,, (M) (M<j/yr) -8-7 . -8 -5.5 -.5 



10 



12 



X 



-A-y ■ -.6 -5.5 -.5 

ni • i i — n~ 



Weak ogfcl/ 
15 fj.m , 
excessc j 





» \, , 

NmAYNLtm 10000 










" ° T^-n 


\ Xi/ysooo 












X^ULoJaoof 12 






v \l(LJ2000 










A x 


v Xlojiooo 






\ L(L o )500 


"m 



10 



(K.-[15]) 



Fig. 9. K s o/(K s -[15])o magnitude-colour diagram of 
ISOGAL sources with DENIS/2MASS counterparts. The 
symbols and the curves are same as shown in Fig. 8. The 
red IRAS sources of Table 3 are denoted by "I" (see text). 
They are also MSX D sources associated with the strong 
15 \im sources. The RGB tip (K s0 ~ 8.0) is shown as dot- 
ted line. Other two regions are same as in Fig. 8. 
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Fig. 10. K s -[15] colour distribution of ISOGAL sources. 
The mass-loss scale at the top of the figure are same as 
shown in Fig. 8. The full line histogram represents the 
stars of all nine fields of our sample (Table 1). The dashed 
line represents stars in the fields with \b\ > 2°. 
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Fig. 11. [15]/[7]-[15] magnitude-colour diagram of 
ISOGAL sources detected in the ISOGAL fields. The 
open circles show the ISOGAL sources having 7/15 [im 
association quality flags of 3 & 4. The crosses show the 
sources having association quality flags of 2. The filled 
hexagons display the associations with 2MASS (since the 
DENIS sources [K s < 6.5] are saturated, they have been 
replaced by 2MASS). The remaining saturated DENIS 
sources (which have no 2MASS counterparts) are shown 
by filled squares. The bright ISOGAL sources associated 
with 2MASS sources (which have no DENIS counterparts) 
are shown by open triangles. The filled triangles show 
the ISOGAL sources which have no available DENIS nor 
2MASS counterparts, mostly because of saturation or 
non observation. The asterisks show the sources detected 
at 8 /im in the MSX survey which are associated with 
strong ISOGAL sources ([15] < 4.5 mag). Among them 
the dark ones are those associated with the MSX D (15 
\xm) sources ([15] < 3.5 mag). The IRAS sources are 
shown by open squares and the brightest ones are denoted 
as "I" (sec text). Sources with larger luminosities and 
super-wind mass-loss characteristic of Miras are inside 
the boundaries shown as dotted lines. Other two regions 
are same as in Fig. 8. 



Fig. 12. [7]o/(K s -[7])o magnitude-colour diagram of 
ISOGAL sources with DENIS/2MASS counterparts. The 
open circles show the sources having ISOGAL-DENIS as- 
sociation quality flags of 4 & 5. The crosses show the 
sources having association quality flags of 2 & 3 (see the 
text). The filled hexagons display the ISOGAL sources 
associated with 2MASS counterparts (since the DENIS 
sources (K s < 6.5) are saturated, they have been replaced 
by the 2MASS sources) . The remaining saturated DENIS 
sources (which have no 2MASS counterparts) are shown 
by filled squares. The bright ISOGAL sources associated 
with 2MASS sources (which have no DENIS counterparts) 
are shown by open triangles. Sources without 15 [ira as- 
sociations are displayed by filled circles. The filled trian- 
gles show the ISOGAL sources associated with 2MASS 
sources but having no DENIS counterparts. The asterisks 
show the sources detected at 8 fim in the MSX survey 
which are associated with strong ISOGAL sources ([15] < 
4.5 mag). Among them the dark ones are those associated 
with the MSX D (15 fim) sources ([15] < 3.5 mag). Stars 
11-4 are outside of the figure limits with (K s -[7])o between 
4 and 10 (see Table 3). The region of high mass- loss AGB 
stars ((K s -[7])o > 1, [7]o < 7) is drawn in the upper right 
corner (see Fig. 9 of Omont et al. 1999). 



to further analysis and observations of mid-IR spectra, an 
essential element for the accurate calibration of the rela- 
tion between M and (K s -[15])o is the detailed theoretical 
analysis of mass-loss physics, dust formation and radia- 
tive transfer, such as performed by Jeong et al. (2002 and 
in preparation). The latter is based on a consistent time 
dependent treatment of hydrodynamics, thermodynamics, 
equilibrium chemistry and dust formation. Such a theoret- 
ical modelling is presently the only way to infer the gas-to- 
dust ratio and the mass-loss rates for such bulge sources 
in the absence of any direct measurement of the mass of 



circumstellar gas and even of the expansion velocity. As 
they quote, the relation derived between M and (K-[12]) 
by Jeong et al. (2002) gives substantially higher mass-loss 
rates than the empirical relation of Lc Bertre and Winters 
(1998) (Fig. 13). On the other hand, the agreement is bet- 
ter with the empirical values of Whitelock et al. (1994) de- 
rived from the classical prescription of Jura (1987) giving 
M from far- infrared emission (Fig. 13). Here, we will apply 
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the explicit relation between M and (K-[15])o derived by 
Jeong ct al. (2002), 



log M = -6.83/(K - [15]) - 3.78 



(1) 



with M in M Q /yr (Fig. 13). 

Alternatively, one could also infer M from (K s -[7])o or 
from [7]-[15]. Deriving M from [7]- [15] (Alard et al. 2001) 
has the advantage that the latter quantity depends very 
little on the extinction. However, this method is less sensi- 
tive since [7]- [15] varies little. On the other hand, (K s -[7])o 
has a low sensitivity for small mass-loss rates, but could 
be conveniently used for large ones. 

It is important to notice that there may also be an 
appreciable uncertainty in the determination of mass-loss 
rates from (K-[15])o for strongly variable stars where we 
have only single epoch measurements of K s and [15], and 
at different epochs for K s and [15] (Table 1). As quoted, 
this uncertainty is reduced by taking the average of DENIS 
and 2MASS K s values and of ISOGAL and MSX 15 [im 
magnitudes. One can thus estimate that the global r.m.s. 
uncertainty on (K s -[15])o remains smaller than 0.3-0.4 
mag, except for the weakest sources (K s0 > 9.0, see Sect. 
6.2). However, because of the rapid variation of M versus 
(K s -[15])o (Fig. 13), the corresponding observational un- 
certainty on M ranges between a factor ~2 for its largest 
values > 10~ 6 M /yr, and ~5 for the smallest ones < 10~ 7 
M©/yr. Such values do not include the model uncertainty 
in the relation between M and (K s -[15])o- 



3 1-6 




Fig. 13. Mass-loss rate (M) as a function of colour. For 
X = (K s -[15])o, Eq. (1) recommended by Jeong et al. 
(2002, JWBS02) is used. For X = (K-[12]) , the three 
curves drawn are derived from Jeong et al. (2002) with 
logio M = -7.69/(K-[12]) ) - 3.40, Whitclock ct al. (1994, 
WMF94) and Le Bertre & Winters (1998, LBW98), re- 
spectively. 



8. Total mass-loss rate in the bulge 

Fig. 14a shows the number distribution of mass-loss rates 
(M) of ISOGAL sources in the bulge as a function of M 
from Eq. (1) along models of Jeong et al. (2002). Fig. 
14b displays the corresponding values of the average total 
mass-loss rate per square degree and per 0.5 bin of log M. 
It should be remembered that the number of sources is 
incomplete for M ^ 10~ 7 M Q /yr. It is not excluded that 
still lower mass- loss rates, M < 10~ 8 M Q /yr, contribute 
a little to the integrated mass-loss in the bulge, because 
of the very large number of stars with 1 (K s -[15])o Ss 2 
(Fig. 10), where such 15 nm excess, if real, still shows the 
presence of circumstellar dust. 

However, because of the very steep distribution of the 
number of sources versus (K s -[15])o (Fig. 10) and the 
rapid variation of M versus (K s -[15])o (Fig. 13), the er- 
ror in (K s -[15])o due to strong variability makes the con- 
tribution of the corresponding bins larger than in reality. 
Therefore, in Figs. 14b & d, we have reduced the contribu- 
tions of the three bins between 10~ 7 and 3 10~ 6 M Q yr _1 
by factors 1.4, 1.2 and 1.1, respectively, as shown in Figs. 
14a & c. However, we have not corrected the bins below 
10~ 7 M Q yr _1 , because the incompleteness of detections of 
such stars could compensate for this effect. 

It is difficult to fully appreciate the uncertainty in the 
data of Figs. 14a & 14b. There are various sources of ob- 
servational uncertainty: errors in K s -[15])o and their un- 
certain correction, bulge membership, Poisson error due 



Table 4. Integrated mass-loss (M Q /yr/deg 2 ) in the "in- 
termediate" bulge, from Eq. (1) (Jeong et al. 2002) with 
M > 10~ 8 M /yr 





with I1-4 T 


without I1-4 1 " 


All fields 


8.1xl(T 4 


3.7xl(T 4 


Fields with \b\ > 2° 


1.2xl(T 4 


1.2xl(T 4 


Fields with \b\ < 1.5° 


1.4xl(T 3 


5.9xl(T 4 



Stars with M 10 5 M©/yr whose bulge membership is 
uncertain (Table 3 & Sect. 6.1) 



to the small number of sources with large M, incomplete- 
ness of detections of stars with small M, etc. They could 
combine to give a substantial error, but probably smaller 
than a factor two. However, the major source of error 
is probably still the calibration of the relation between 
K s -[15])o & M, which depends on dust properties and the 
dust-to-gas ratio. Its modelling depends on complex stel- 
lar properties and remains difficult despite the spectacular 
progress achieved by Jeong et al. (2002). 

It is seen in Fig. 14b that the contribution to the inte- 
grated mass-loss per bin of log M is dominated by mass- 
loss rates larger than 10 -6 M Q /yr. It is even possible that 
the four stars 11-4 with (K s -[15]) > 6, M ^ 10~ 5 M Q /yr 
(Table 3), contribute for about half the total. But their 
belonging to the bulge is not fully established (see Sect. 
6). 
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Numerical values for the integrated mass-loss are dis- 
played in Table 4. First are values for mass-loss per 
square degree averaged over all nine fields of Table 1, 
with and without including the contribution of the four 
stars 11-4, respectively. Since it is difficult to decide 
whether such stars belong to the central bulge/disk or 
not (Sect. 6.1), we also give the same quantities for the 
three fields with |6| > 2° where such stars with very large 
values of M ((K s -[15])o > 4) are absent (see Fig. 10). 
The corresponding distribution of mass-loss rates are dis- 
played in Figs. 14c and 14d. In the absence of stars with 
(K s -[15])o > 4, it is seen that the relative contribution 
of large mass-loss rates to the integrated mass-loss seems 
smaller. However, the total number of relevant stars with 
(K s -[15])o > 2 is about four times smaller in the three 
fields with |6| > 2° than in the totality of the fields. 
Therefore, while there are only four stars with M > 1CP 5 
M©/yr in the totality of the fields, one cannot consider 
that their absence in fields with \b\ > 2° is statistically 
significant. 

Such results seems to show that the relative contribu- 
tion of stars with low mass-loss rates (M < 1CT 6 M Q /yr) 
is more important in the mass return to the interstel- 
lar medium for the old bulge population, which is more 
predominant at \b\ > 2°. Even the stars with M < 10~ 7 
M /yr may have a significant contribution. 

It is not clear whether the stars with very large mass- 
loss rates, M i> 10 -5 M Q /yr, which are probably younger 
with a relatively large initial mass, may have an essential 
contribution, as it seems the case for the solar neighbour- 
hood (see e.g. Le Bertre et al. 2001). It would be interest- 
ing to make a detailed comparison with the solar neigh- 
bourhood. However, in addition to the question of bulge 
membership of stars with such large mass-loss in our sam- 
ple, it is difficult to make a precise comparison with the 
recent work of Le Bertre et al. (2001) because their sample 
is not as complete for stars with low mass-loss as for large 
mass-loss ones. 

We may assume the value of Eq. (A. 8) for an approxi- 
mate average bulge projected mass density at the average 
position I = 0°, |6| = 2° for our fields, from the model of 
Wainscoat et al. (1992) - £ Bu (0,2) ~ 5.7xl0 7 M© dcg~ 2 . 
We infer a rough value for the total (bulge + disk) pro- 
jected mass density of ~ 8xl0 7 M© deg~ 2 , by assuming 
the same disk/bulge ratio, 30/70, as for the stellar counts 
estimated in the conclusion of Sect. 4.2. The derived inte- 
grated mass-loss rate per unit stellar mass is thus about 
0.4-1.0 10" 11 yr- 1 . 

One notes that there is a difference by a large factor in 
the average mass-loss per square degree between the fields 
with \b\ > 2° and those with |6| < 1.5°. This factor is 5 
and 12, without and with the contribution of stars 11-4, 
respectively (Table 4). On the other hand, since the bulge 
projected mass density varies as [b 2 +(£/1.6) 2 ]~ 0A from the 
model of Wainscoat et al. (1992) (see Eq. A. 7)), the cor- 
responding ratio for the average projected mass density is 
only ~ 2 for the bulge. It should be ~ 3 when one takes 
into account the disk contribution. A possible explanation 



is that the fields closer to the Galactic Center contain a 
larger proportion of younger stars than in the outer bulge 
(Frogcl et al. 1999), which, in particular, contribute more 
to the large mass-loss rates. However, a deeper discussion, 
with more elaborate Galactic models and a better iden- 
tification of foreground sources, will be necessary to well 
establish and quantify these properties. 

9. Conclusion 

We have analysed in detail the AGB population in the 
available ISOGAL "intermediate" bulge fields (|Z| < 2°, 
|6| ~ 1° - 4°). We have confirmed that the combination of 
near-infrared and mid-infrared (7 and 15 fim) ISOGAL 
data allows reliable detection of AGB stars above the 
RGB tip with a determination of their luminosity (pro- 
viding that they belong to the bulge) and mass-loss rate. 
We conclude that almost all the ^1300 ISOGAL sources 
detected both at 7 and 15 firn on the line of sight of 
the bulge are AGB stars or RGB tip stars. A large pro- 
portion of these AGB stars have appreciable mass-loss 
rates, as shown by the excess in (K s -[15])o and ([7]-[15])o 
colours, characteristic of circumstcllar dust emission. We 
have performed a preliminary determination of mass-loss 
rates from (K s -[15])o from recent theoretical modelling by 
Jcong et al. (2002) which are in agreement with Whitelock 
et al. (1994) for (K-[12]) . For the bulk of our sample, and 
especially at |6| <~ 1°, the total mass returned to the in- 
terstellar medium is dominated by mass-loss rates larger 
than 10~ 6 M©/yr. However, four stars, out of ~2000, have 
mass-loss rates 10~ 5 M©/yr. They could dominate the 
mass return to the interstellar medium if they all belonged 
to the bulge, which is still unclear. In the more outer bulge, 
|6| > 2°, the contribution of large mass-loss rates seems 
less important. 

Appendix A: Numbers of foreground and 
background disk sources 

In order to derive an order of magnitude of the propor- 
tion of foreground and background disk stars, on the line 
of sight of the observed bulge fields, one can use any rea- 
sonable simple Galactic model for the distribution of con- 
tributing stars, e.g. the one by Wainscoat et al. (1992). We 
have quoted that practically all bulge sources detected by 
ISOGAL are M giants later than about M3. Background 
detected stars should be even later M giants, while the 
great majority of foreground disk sources are K or early 
M giants. It would be straightforward, but tedious, to de- 
termine the proportion of foreground and background disk 
stars in any given magnitude bin, on various lines of sight, 
by direct numerical integrations of the Galactic distribu- 
tions of the various classes of stars - especially red giants 
- of Wainscoat et al. However, we take advantage of the 
simple behaviour of the giant luminosity function to make 
an approximate analytical integration. The K-M giant lu- 
minosity function is roughly proportional to L _1 (see e.g. 
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Fig. 14. (a). Distribution of mass-loss rates (M) of ISOGAL sources inferred from (K s -[15])o in the nine bulge holds, 
derived from the formula of Jeong et al. (2002), Eq. (1). The dotted line includes a tentative correction for the spurious 
counts due to the error in (K s -[15])o and the very steep count distribution, as described in the text. (b). Average total 
mass-loss rate per square degree and per 0.5 bin of log (M) as a function of M after this correction. The contribution 
of stars 11-4, M ^ 10~ 5 M /yr, is shown by a cross. The arrow shows the incompleteness on the data for the smallest 
mass-loss rates, (c) & (d). Same as a & b, but limited to the three fields with |6| > 2. 



Frogel & Whitford 1987). At the point (r,z) in the direc- 
tion (£, b), with b « 1 and z = rb, the number of such 
giants per unit volume and absolute magnitude bin dMi 
is p g (r,z) x ^(M^dMi where fi(M 4 ) is approximately pro- 
portional to 10 M '/ 2 ' 5 . The number of stars per unit volume 
at point (r,z) contributing to the observed magnitude bin 
dmi, with Mi = Mi - 2.51og(r/10pc) 2 , is thus proportional 
to r~ 2 x p s (r,z). The number of counts per unit of solid 



angle ft (e.g. deg 2 ) in a magnitude bin in the direction (£, 
b), is then proportional to 

jr-i^^drotj^^dr (A.l) 

We have therefore integrated the density models of 
Eqs. (4) and (5) of Wainscoat et al. (1992) for the disk 
and the bulge, respectively, on a typical line of sight for 
our fields (e.g. I = 0, b = 2°). For the disk, 
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p gD (r,z) = p Qg exp[-(R- Ro)/h - \z\/h z ] (A.2) 

where p$ g is the disk density of red giants in the solar 
neigbourhood, R is the distance from the Galactic Center 
within the plane, Ro = 8.5 kpc is the distance from the Sun 
to the Galactic Center, h = 3.5 kpc and h z = 325 pc for the 
giants. For r < Ro, R-Ro = -r; r > Ro, R-Ro = r-2Ro; 
and at b = 2°, z/h z = r/(9.3kpc). The integrations of 
Eqs. (A. 1-2) over r along the line of sight (£, b), from to 
Ro and from Ro to oo, respectively, are straightforward, 
yielding, in the direction £—0, 6=2°, Np oc 20p Og for the 
disk counts with r < R , and N B k oc 12p 0g for the disk 
counts with r > R . 

For the bulge, 

p gB {r, z) = 3.6 x p 0g x- 1 - 8 cxp (-x 3 ) (A.3) 
where 



x = y/R 2 + {l.6z) 2 /R{ 

= ^R 2 [£ 2 + (1.66) 2 ] + r' 2 /R 2 (A.4) 

where Ri = 2.0 kpc is the bulge "radius". However, 
for the small values of £ & b in the considered fields, most 
of the integral of p g s(r,z) comes from bulge points with 
x << 1. Therefore, the term exp (— x 3 ) may be neglected 
in Eq. (A.3). In the direction £=0, 6=2°, the bulge counts 
are thus proportional to 

/+oo 
x~ 1& dr' (A.5) 
-OO 

or 

3.6 x POg R^Ro/R^lb 2 + (£/1.6) 2 ]- 0A 

/+oo 
(y 2 + l)- - 9 dy (A.6) 
-oo 

yielding for the bulge counts, in the direction £=0, 
6=2°, N Bu ex 57p 0g . 

With this model, on the line of sight £=0, 6=2°, one 
thus expects to find about 65% of bulge sources, 22% of 
disk sources with r < R , and 13% of disk sources with 
r > R , among the ISOGAL detected sources. 

Remark. Similar integrals intervene in the com- 
putation of the total stellar mass column densities 
M c d = Jq°° Pg{ r , z ) dr, yielding, with this model: 

M cdBu (M Q kpc- 2 ) = 

100 x p M0 (MQkpc- 3 ) x [b 2 + (£/1.6) 2 }~ 0A (A.7) 



for the mass column density of the bulge, where b 
and £ are expressed in degrees, puo is the stellar mass 
density in the solar neighbourhood expressed in Mq 
per kpc 3 . We will adopt the value of Chabricr (2001) 
Pmo = 4.5 10 7 M Q fcpc~ 3 . Since most of the bulge sources 
are close to the Galactic Center, one can directly infer the 
projected bulge mass per solid angle unit, e.g. per square 
degree. For instance, in the direction £=0, b~2°, 

S B u(o,2) « 5.7 x 10 7 M Q deg- 2 (A.8) 
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